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Graphene is a realization of an esoteric class of materials - electronic crystalline membranes. We 
study the interplay between the free electrons and the two-dimensional crystal, and find that it in- 
duces a substantial effect on the elastic structure of the membrane. For the hole-doped membrane, 
in particular, we predict a spontaneous buckling. In addition, attenuation of elastic waves is ex- 
pected, due to the effect of corrugations on the bulk modulus. These discoveries have a considerable 
magnitude in graphene, affecting both its mesoscopic structure, and its electrical resistivity, which 
has an inherent asymmetry between hole- and electron-doped graphene. 
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I. INTRODUCTION 



II. ELASTIC PROPERTIES OF GRAPHENE 



The isolation of graphene^ has no t only provided the 
catalyst for countless discoverie^^is]^ ^^y^^ ^j^gQ heralded 
the unification of two independent branches of phys ics: 
the structure of crystalline membranes on one hand,!^ 
and the properties of quantum field theories (QFTs) of 
fermions in two dimensions (2D) on the other hand.'^ In 
its structure, graphene is the ultimate crystalline mem- 
brane - made of carbon atoms arranged in a 2D hexag- 
onal lattice. Its low-energy electronic structure is analo- 
gous to a massless Dirac fermion, establishing graphene 
as a table-top experimental device for the study of 2D 
QFT.'S^ Clearly, graphene offers an unprecedented oppor- 
tunity to examine the interplay between these aspects of 
the material. 

The two-dimensional structure is known to be an impor- 
tant element in the electronic transport properties, due 
to scattering of ch arge c arriers off corrugation^ and in- 
plane deformation^i^^. However, is the surface struc- 
ture affected by the electronic properties? In this Brief 
Report we argue that the answer to this question is posi- 
tive. We evaluate the contribution of the free-electrons to 
the elastic free energy, and find that the electron-phonon 
interactions lead to a substantial effect on the elastic 
constants. In particular, we find that the bulk mod- 
ulus acquires a correction which depends on the wave- 
length of the excited elastic waves, thus implying an at- 
tenuation of these waves. In addition, and more impor- 
tant, we find that a mesoscopic appearance of ripples, 
i.e., buckling, is possible for a hole-doped membrane. To 
our knowledge, this is a unique property of electronic 
crystalline membranes in general, and graphene in par- 
ticular. Seeing that the scattering of charge-carriers off 
ripples increases the resistivity of the material, this pre- 
diction hints that the origin of the asymmetry found 
experimentalljli^lisi between the resistivity of hole- and 
electron-doped graphene, is this coupling of the physical 
structure and electronic characteristics. As such, this es- 
oteric effect should be taken into account when designing 
graphene based technological applications. 



As a 2D material, the mere existence of graphene 
was in doubt, due to the Mermin- Wagner theorem,!^ 
which forbids the existence of long-range order in two- 
dimensions. This seeming contradiction can be resolved 
by introducing small out- of-pla ne crumpling, which sup- 
presses thermal vibrations .SUSl Extensive theoretical, nu- 
merical, and experimental investigations of crystalline 
membranes have established the existence and stability of 
a low temperature phase, characterized b^ corrugations 
and ripples, however asymptotically flat.^ The stability 
of graphene, even without the support of a substrate, was 
found experimentalljff^ and numerically.^''^ 

In order to model this "almost-flat" phase of a mem- 
brane, it is useful to describe the deviation from the 
ideal phase, i.e., a flat surface with perfect lattice. In- 
plane deformations are characterized by a two dimen- 
sional vector fleld u, and the out-of-plane deformation by 
a fleld h. This deformation changes the distance between 
two points on the surface, initially separated by an in- 
flnitesimal vector dx — {dxi , dx2), by 2 ■ 
where Uij is the strain tensor, Uij — }^{diUj + djUi) + 
\{dih){djh). The mesoscopic structure of the surface is 
determined by a free energy functional, which has to pre- 
serve the homogeneity and isotropy of this limit - the 
elastic free energy: 



F[u, h] 



d^XK{Ahf + (1) 



d^x 



i,j = l 



The coefficients are the elastic properties of the mem- 
brane, its bending energy k « 1.1 eV, bulk modulus 
« 11 eV A~^, and shear modulus « 9 eV A~^]i^^ 
They originate in the a bond between the carbon atoms, 
which is a consequence of sp^ hybridization, that forms 
a deep valence band. This fact is evident in the high cre- 
ation energy for structural defects, reflected experimen- 



2 



tally in the high lattice quality. In fact, lattice defects 
were not observe d in graphene, even at strain values of 
about validating the elastic approximation up 

to high strains. 



III. ELECTRONIC STRUCTURE OF 
CORRUGATED GRAPHENE 

An additional electron in a p orbital differentiates 
graphene from other crystalline membranes. This orbital, 
which is perpendicular to the planar structure, forms a 
half filled tt band. This electron is responsible for one 
of the unusual aspects of graphene ~ its celebrated Dirac 
type low-energy excitations. At low-energy, the hexago- 
nal lattice leads to a linearly vanishing density of states 
around two points, named Dirac points, in the recipro- 
cal space This dispersion relation characterizes mass- 
less, chiral Dirac fermions, with an effective speed of light 
Vf ^ c/300, whose size is dictated by the hopping inte- 
gral of the carbon bond. This picture leads to a ballistic 
movement of the electron. Experimentally this mobil- 
ity is limited, probably due to the interaction of electro n 
with phonon fields, i.e., the surface structure .'^^'^3121111 

The dominant effect the complex mesoscopic structure 
of graphene has on the ideal QED electronic structure 
arises from the deformation of the lattice. The source of 
this deformation potential is the local change in the Fermi 
energy measured from the bottom of the a bond, propor- 
tional to the change in area SS cx X]?=i ''^a- The result- 
ing deformation potential has the form Vg = ^(X]i=i ''^a) 
with, D = 10 — 30 eV.'^We note, as stated by Hwang and 
Das-Sarma,^° that since the source of this potential is not 
a Coulomb potential, its screening can be neglected. 

Another effect of the deformation of the lattice is a 
change in the hoping integral, which ma nifests itself as 
an effective gauge field in the Dirac picture.'^i^Its struc- 
ture is determined by the symmetry of the lattice, and 
its characteristic size is much smaller than the deforma- 
tion energy, 52 ~ 1 — 4eV,'2ilthus we will postpone the 
discussion on its effect on the elastic free enegy to a dif- 
ferent work.[22l We comment, however, that this pesudo 
magnetic field leads to preferred directions of ripples, fol- 
lowing the underlying lattice symmetry. 

The energetics of the tt electron, in the continuum 
limit, can be summarized in a QED like action, of a 
fermion in non-zero chemical potential. The Matsubara 
action of an electron, of inverse temperature /3, is 



electro-chemical potential V is a sum of the deformation 
energy, and of any other external electro-chemical po- 
tential SV. As will be shown later, the interaction of the 
deformation energy with an external electro-chemical po- 
tential can lead to buckling in graphene, as the external 
potential can compensate a deformation energy, allow- 
ing a deformed lattice. In this way, a buckled phase can 
be favored. In order for this state to be the equilibrium 
phase, a mutual minimization of the elastic me mbran e 
energy (Eq. [l]) , and the electron energy is needed.l^^iail 



IV. CONTRIBUTION OF FREE ELECTRONS 
TO THE FREE ENERGY 

Prominent is the fact that the characteristic time scale 
of the electrons is determined by Vf sa 10^ cm/sec, 
whereas the mesoscopic structure of graphene changes 
in a much longer time scale, dictated by the speed of 
sound Vph « 2 • 10^ cm/sec. As a result, electronic excita- 
tions are substantially more energetic than the phonons, 
or elastic, excitations, and thus the electronic degrees- 
of-freedom can be integrated out. In addition, we as- 
sume that the strength of the electro-chemical potential is 
weak, keeping terms only to second order in the strength 
V^. The resulting contribution of the tt electrons to the 

free energy is Feiectron = \j (^n^(g; I3)V^ , where the 
subscript q denotes Fourier transform, i.e., for a constant 
external electro-chemical potential: Vq = £^(X)i=i + 
{2ttY5V 5{q). n^(g; (3) is the polarization operator of mo- 
mentum g, which in the one-loop level is: 



2 n— 00 



71 — — 00 



d?k ka{k + q)p ^ ^ 
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where the Matsubara frequencies are k^ — ^("' + 5)- 
The polarization operator of the static electro-chemical 
potential, i.e., q'^ 



n^(q;/3) 



0, is thus.*^ 

dxln 



2cosh(^y.?(r3^ 



7r(fiu/)2/3 Jo 

For hvfPq 1, one can show that n^(g;/3) — ^^-^^"^2^ . 
In the opposite limit hvf[3q 3> 1, n^(g;/3) = -^j^ 
The transition between these asymptotic behaviors oc- 
curs rather abruptly, for hvf(3q ~ 2, as can be seen in 
Figure [T] 



■ {j°ido + t{Duu + SV)) + vn-V^ •«'n(T,f). 

Here, n is the index due to the N — 2 spin degeneracy, 
^ is a four component Dirac fermion, reflecting the two 
sub-lattices and two Dirac cones. The 7 matrices satisfy 
the Clifford algebra, {7^, jv} = 25 Note that the total 



V. EFFECT OF FREE ELECTRONS ON THE 
BULK MODULUS 

One can clearly see that the contribution to the free 
energy leads to a term in the free energy which is pro- 
portional to I X]i=i ''^iiP' i-^-; increases the bulk modu- 
lus. Specifically, the Lame constant A is renormalized, 
\{q) = A -f D^I\^ {q\ (3). This implies a finite contribu- 
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FIG. 1: (Color online) Dimensionless polarization operator 
{hvf)^PIl^ as a function of the dimensionless momentum 
hvfPq. 



ripples on the order o f 50 — 150 A, which are found experi- 
mentally on graphen j^^ l ^^ ^and around room temperature, 
the appropriate limit is hvfPq ^ 1, which implies an 

addition of 3.85 eV • f^)^ {~~A^ 



This 



\30cVJ 

tribution is substantial compared to the experimentally 
found X — 2 eV A , and depends on the wave-number, 
thus should be observed as an attenuation of sound wave 
packets in graphene. 

The dependence of this attenuation on the deforma- 
tion energy D is quadratic. This strong dependence can 
be used to measure this parameter, whose value is known 
rather poorly. This parameter is of immense importance, 
as the resistivity induced by the deformation potential 
is considered to be the dominant source of carrier scat- 
tering at high carrier densities.^ D is usually extracted 
from its contribution to resistivity, hence an independent 
measurement, as suggested here, is valuable. 



VI. SPONTANEOUS BUCKLING 

However, the most interesting effect stems from the 
term in F^iectron proportional to the cross-product of 5Vq 

and ^iij which has the form 6F = t J d^aJ'X^i^i "^a 

with r = DU^{q = 0; f3)SV. This term leads to a buck- 
ling instability. The reason for that is easily seen when 
considering vanishing in-plane deformations, and keeping 
quadratic terms in h. In this case: 



ions, which are held by the elastic forces. The competi- 
tion can be resolved by the ripple creation. A positive 
potential, however, results in local stretching, with no 
advantage in rippling. 

The external potential SV, in the single electron pic- 
ture, is reflected in the chemical potential of an electron. 
Sources of such chemical potential are numerous. For 
example, adsorption of molecules on graphene could pro- 
vide a simple possible mechanism to check the theory, 
predicting an inherent difference between exothermic and 
endothermic adsorptions. For example, doping with NO2 
molecules, which behave as holes (acceptors), has been 
demonstrated to high doping levels, n oc 10^^cm~^ !^^ l ^^ l 
Another interesting source of chemical potential is a fi- 
nite density of surface charge carriers. A source of such 
charge carriers, of immense importance for future tech- 
nological applications, is an external gate, through the 
electric field effect^ The charge carrier density n is a re- 
sult of the capacitance of the gate, thus n is linear in the 
gate voltage Vg. Both in chemical doping and in electri- 
cal doping, the resulting change in the chemical potential 
can be estimated by SV = sgn{n)hv j n\n\ . 

As a result, the buckling is characterized by ripples 



with a characteristic size, ^ = 2tt 
cally: 



4ln2DSV ' 



or numeri- 



^ w I44A 



300°K 



D 



30 eV 



1012cm- 



(3) 

In the case of adsorbed molecules, where carrier concen- 
trations on the order of lOi^cm"^ have been reached,!^ 
we predict ripples of wavelength as small as 100 A. For 
electrical doping, where the carrier concentration is con- 
trolled by the gate voltage, the size of the induced ripples 

is easily calculated to be ^ w 270 A (uw) 

This rather amazing outcome, that can and should be 
checked experimentally, reflects the possibility of geomet- 
ric response of the graphene to external potential. It 
is clear, that this effect is valid mainly for suspended 
graphene, as the forces exerted by a substrate would pin 
the structure of the surface to that of the substrate. In 
addition, it is important that the graphene would not be 
covered by a thin layer of water, as suggested in Ref. 27, 
as this will screen the effect of adsorbants. 



F[u, h] 



k(A/i)2-Ht(V/i)^ 



(2) 



If T < 0, this free-energy is minimized when the spa- 
tial configuration has a correlation length of the order 
of ^ = 27ry^^j^.l^ 22 | 26| Intuitively, one may understand the 

condition r < 0, by thinking about a piece of paper, 
which is not buckled when stretched, but is buckled when 
contracted. The same happens here, with the chemical 
potential, which takes the part of the contraction stress 
in the analogy to a paper. A negative potential pulls the 



VII. ADDED RESISTIVITY INDUCED BY THE 
BUCKLING 

These ripples, predicted for hole-doped gra phene only, 
will increase the resistivity of the matter,'^'!^' creating an 
asymmetry in the transport properties of hole and elec- 
tron doped graphene. Such an asymmetry is also ex- 
pected due to charged impurities in the substrate ,1221 and 
due to a formation of a p-n junction at the contacts.!^ 
In order to estimate the contribution of corrugations to 
the electrical resistivity, we use the calculations of Ref. 9, 
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where this added resistivity is related to the scattering 
of electrons on excited flexural phonons within the rip- 
ples. For low temperatures k^T < 2tt^^ k. 100°K, 
phonon excitations are not allowed, thus these will not 
contribute to the resistivity. For high temperatures, one 
can approximate the contribution to the resistivity by 

^ 200 - 600f^.(^).(J^) 

where a is the interatomic distance. 

As previously mentioned, corrugations on the order of 
50 - 150 A were found in suspended graphene sheets.'^ 
This fact can be explained by the effect presented here, 
through molecular hole-doping. However, the hole-doped 
buckling can not reproduce the buckling found in quan- 
tum Monte-Carlo simulations, which do not take into ac- 
count the TT-electrons.l^In Refs. 16 and 17, these ripples 
were associated with quenched-ripples originating in a 
substrate, or with peculiarities in the carbon bond. The 
production of ripples in the process of isolation of one 
graphene layer, which competes with the effect presented 
here, can be minimized by depositing graphene on liq- 
uid substrates. ^"^ However, buckling which originates in 
the carbon bond will manifest itself, in the elastic free 
energy, with the same term that leads to the effect dis- 
cussed here.'^Sl As a result it will interfere constructively 
with the effect predicted here. A distinction between all 
these effects can be made by the different dependence on 
the thermodynamic conditions, as the contribution to re- 
sistivity of quenched ripples or those originating in the 
properties of the chemical bond are quadratic in tem- 
perature and inverse proportional to the charge carrier 
concentratiorF^ and do not create an electron- hole asym- 
metry. 

The asymmetry in the electron mobility between hole 
doped and electron doped samples has been found in 
many experiments, e.g. Refs. 12 and 13. In particu- 
lar, recently, Bolotin et aP^ pointed out an unexpected 
dependence of the resistivity in charge carrier densities, 
mainly in hole-doped graphene, that can be qualitatively 
explained by the effect predicted in the current work.l^ 



VIII. SUMMARY 

We have presented a theoretical investigation of the 
effect 2D relativistic fermions have on the structure of 
the 2D lattice on which they reside. This is achieved by 
evaluating the contribution of the fermions to the elastic 
free energy of the surface, a method which leads to the 
dependence of the elastic constants on the momentum 
of the excitations. This phenomenon is interpreted as 
an attenuation of sound waves excited on the surface. 
In addition, we predict a possibility of structural change 
reflected in a spontaneous buckling of the surface, due to 
a non-zero fermion chemical potential. 

This study is motivated by, and demonstrated on, 
graphene. It would be interesting to see whether the 
same theoretical approach can lead to the better under- 
standing of other materials of two-dimensional character. 
In graphene, these phenomena are found to have an ob- 
servable effect on the elastic constants, due to the large 
deformation energy of the lattice, thus offering a different 
experimental path to measure this poorly-known quan- 
tity. The spontaneous buckling of hole doped graphene, 
which can be realized by chemical doping or by electri- 
cal doping via the electric-field effect, is expected to cre- 
ate ripples which can be directly observed, and to have 
an effect on the resistivity of this phase of the material. 
Clearly, an experiment which will provide a direct ob- 
servation of the rippling as well as the predicted sound 
attenuation, is called for. 

The theoretical procedure presented here provides 
an important tool for characterizing the structure of 
graphene and its transport properties - the key elements 
in the design and quality control of any future technolog- 
ical application. 
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